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ABSTRACT 

The role of binary progenitors of neutron stars in the apparent distribution 
of space velocities and spin-velocity alignment observed in young pulsars is 
studied. A Monte-Carlo synthesis of pulsar population from single and binary 
stars with different assumptions about the NS natal kick model (direction 
distribution, amplitude, and kick reduction in binary progenitors which ex- 
perienced mass exchange due to Roche lobe overflow with initial masses on 
the main sequence from the range 8-11 Mq) is performed. The calculated 
spin-velocity alignment distributions are compared with observational data 
obtained from radio polarization measurements. The observed space velocity 
of pulsars is found to be mostly shaped by the natal kick velocity form and 
its amplitude; the fraction of binaries is not important here for reasonably 
large kicks. The distribution of kick direction relative to the spin axis during 
the formation of a NS is found to affect strongly the spin-velocity correlation 
of pulsars. Comparison with observed pulsar spin- velocity angles favours a 
sizeable fraction of binary progenitors and the kick-spin angle ~ 5 — 20°. The 
form of the initial binary mass ratio distribution does not affect our results. 
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1 INTRODUCTION 



Neutron stars (NS) resulting from core collapses of massive stars are observed in very dif- 
ferent astrophysical sources (radio pulsars, compact central objects in supernova remnants, 
soft gamma-ray repeaters, anomalous X-ray pulsars, cooling radio-quiet NSs, X-ray binaries, 
etc.), and new observational appearences of NSs are discovered nearly every few years (ro- 
tating radio transients are one of the latest examples; f or a comprehensiv e review on recent 



Bassa et al. 



20081 ). Determination of 



studies in this field see, e.g., the conference proceedings 
the birth properties of NSs is crucial to understand the physics of supernovae and to model 
the evolution of NSs. For a given equation of state of the NS matter, main physical param- 
eters of a newborn NS which determine its evolution and observational appearences include 
its mass, spin period, surface magnetic field, and the initial space velocity. The situation is 
more complicated when a NS is formed in a binary system. Then the subsequent evolution 
of the NS will be also determined by the binary orbital parameters and the evolutionary 
state of the secondary component prior to the supernova explosion, in addition to initial 
conditions of the collapse. 

In single radio pulsars, NS space velocities (more precisely, their tangent component) 
can be measured directly. Using complicated models, one can reconstru c t the distribution 
of the NS birth velocities (see, for example, IFaucher-Giguere fc Kaspil (120061 ) and refer- 
ences therein). The shape of the result ing distribution is still ambiguous. Some authors 



( Arzoumanian. Chernoff fc Cordes 



(IHobbs et al. 



2002) favour a bimodal Maxwellian distribution, some 



2005 ) argue that the data can be explained with a single Maxwellian curve, 



while others (IFaucher-Giguere fc Kaspil 120061 ) prefer single-mode non-Maxwellian distribu- 



tions. Moreover, the space velocity (both in magnitude and direction) can be correlated 
with other parameters. In pulsars, one of the best known correlations is found between the 
direction of the space velocity vector and the spin axis of the NS. It is well established for 
the Crab and Vela pulsars and is f ound in several other objects with well studied pulsar 



wind nebulae (jNg fc Romani 



2004 



20071 1. Polarimetric measurements of radio pulses also 



allow to estimate the angle 5 betw een the spin axis and the NS space velocity direction 



(Uohnston et al 



2005 



Rankin 



20071 ). The spin- velocity correlation is firmly established only 



in young pulsars, since for older ones the motion in the galactic gravitational potential 
smears out any initial relation between the NS spin and velocity vectors. 

Recently, 



Johnston et al. 



(120071 ) reported new measurements of the angle 5 between the 
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spin axis and space velocity direction in 14 radio pulsars. Half of them demonstrate strong 
alignment, while others do not. New pulsar measurements and observ ations of puls ar wind 



Ranki 



2003). 



nebulae are expected to increase the list of objects with known 5 (e.g. 

The origin of the spin-velocity alignment is usually related to the kick mechanism oper- 
ating during or immediately after the stellar core collapse. The reason why NSs obtain so 
large (up to > 1000 km s _1 ) space velocities remains unknown. Se yeral absolutely dif ferent 



physical mech anisms are discussed (see, for exa mple, the table in iBombaci fe Popov 



and reviews in 



Lai. Chernoff fc Cordes 



2001 



and 



Postnov fc Yungelson 



2004 



200J). 



Historically, the first proposed kick mechanism was r elated to an asymmetric mass ejec- 
tion during the supernova explosion (jShklovskiil 119701). Recently, a modificati on of this 



Schecketal 



(12004 . 



20061 ). In these 



mechanism was studied in numerical simulations by 
simulations large kicks can be recovered, but the direction of the velocity is found to be 
poorly correlated with the NS spin axis. Nevertheless, in the framework of this mecha- 
nism the correlation is possible for small spin periods of newborn NSs (P < 1 msec), as 
rapid rotation results in avera ging of momentum along all directions except the spin vector 



(ILai. Chernoff fc Cordes 



200 1) , and anisotropics appear not randomly but preferably along 



the rotation axis (jScheck et al.l 120061 ) . 

In some models, the NS kick is related to an asymmetric n eutrino e missio n in a strong 
magnetic field. The first scenario of this type was proposed by IChugail (119841 ) (see also the 



pa per by 



Bisnovatvi-Kogan 



1993). Several variants of this mechanism were later discussed 



by lLai. Chernoff fc Cordes! (120011 1 . In the neutrino scenarios, the spin- velocity alignment is 



possible if the neutrino emission time is longer than the spin period or if the magnetic dipole 
is initially aligned with the spin axis. 

A tight spin- velocity alignment naturally emerges in two other mechanisms: when the kick 



is d ue to asymmetry between two opposite 



NS (IKhokhlov et al. 



1999 



Akiyama et al. 



dire cted matter jets emanating from a newborn 



20031 ). or when a rotating NS with asymmetric 



(e.g. off-centered) magnetic dipole is gradually ac celerated by electromagneti c radiation 



Harrison & Tademaru 



19751 ) . The spin- 



(the electromagnetic rocket mechanism suggested by 
yelocity alignment may be a consequence of th e magneto-rotational supernova explosion 
(lArdeljan. Bisnovatyi-Kogan fc Moiseenkd 120051 1 . 

A perfect orthogonal spin-velocity misalignm ent is expected in the case of the col- 



lapsing core fragmentation due to rapid rotation (IBerezinsky et al 



Imshennik fc Nadezhin 



1992 



1988 : 



Imshennik 



1992 : 



Colpi fc Wassermanll2002l ). However, it is unclear whether this 
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mechanism can be generic for the majority of NSs, since numerical mo dels of stellar evo 



suggest rather slow rotation of the stellar core prior to the collapse (IHeger et al.ll2005l ). 



ution 



If all NSs were originated from isolated progenitors, deviations from the spin-velocity 
alignment could be related to differences in the kick mechanism or in the NS behavior 
immediately after the formation (the typical space velocity of massive stars - NS-progenitors, 
~ 10 — 30 km s -1 , is much smaller than the NS kick velocities and so can not significantly 
affect the relative orientation of the sp in and velocity axis, see the results below). This 
posibility was studied, for example, by iNg fc Romanil (120071 ). However, at least a half of 
isolated NSs are expected to be born in binaries so the net space velocity acquired by a NS 
formed in a binary system should be the sum of the NS natal kick and the orbital velocity 
of the NS progenitor prior to the binary disruption (plus the velocity of the binary system 
barycenter). Even if the kick mechanism generates a single-mode velocity distribution, the 
additional contribution from progenitor's orbital motion in the binary system can change 
its shape. While the effects of the binary progenitor on the observed NS space velocity 
distribution can be noticeable only for small kicks (of order or smaller than the orbital 
velocity of the collapsing star, typically < 100 km s -1 ), they can appear more pronoun ced 
in the NS spin- velocity misalignment. In our previous paper (IPostnov &: Kuranovl 120081 ) we 
studied the NS spin-kick velocity correlation effect on binary NS coalescence rates and spin- 
orbit misalignment of the components. In the present paper we specially address the role of 
binaries in shaping the NS space velocity distribution and their effect on the spin-velocity 
alignment in observed radio pulsars. 
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2 POPULATION SYNTHESIS CODE 

For our population synthesis calculations we use the SCENARIO MACHINE code, de - 
veloped at Sternberg Astronomical Institute (see iLipunov. Postnov fc Prokhorovl (119961 ); 



Lipunov et al 



(120071 ) and references therein^). The gene r al de scription of the population 
synthesis techniques can be found in iPopov fc Prokhorovl (120071 ) . 

The standard assumptions about the binary evolution are made: the Salpeter mass func- 
tion for the primary mass, dN/dM 1 ~ M-f 2 ' 35 , the initial semi-major axis distribution in 
Opik's form dN/dloga = const, 



10ik 



max 



< a < 10 7 R P 



©) 



R L (M ly 

where R^(Mi) is the radius of the Roche lobe of the primary component. 

Calculations are done for two assumptions about the initial mass ratio (q = M^jMx < 1) 
distribution: dN/dq ~ q° (flat distribution), and the mass ratio distribution dN/dq ~ q 2 
which is strongly peaked to wards equal masses of t he co mponents. The latter choice is 
related to recent studies by iPinsonneault fc Stanekl (120061 ) who found a large fraction of 
twins among 21 eclipsing massive binaries in SMC. The real initial distribution of binary 
mass ratio in our galaxy is still unknown, but the choice of these two limiting cases (the flat 
one and the one strongly peaked towards equal masses) gives us a flavor how sensitive are 
our results to this important initial parameter. 

The common envelope p hase of the binary evolution is treated in the standard way based 
on the energy conservation (IPostnov &: Yungelsonll2006l ) with the efficiency acE = 0.5. 

We assume that the kick velocity vector is confined within a cone coaxial with the 
progenitor's rotation axis and characterized by the angle 9 ^ 7r/2. We shall consider only 



central k icks thus ignoring t 



NS spin ( Spruit fc Phinne 



leoretically feasible off-center kicks simultaneously a 



1998 



Postnov fc Prokhorov 



1998 



Wang. Lai fc Han 



: ecting the 



20071 ). 



We made two assumptions about the NS kicks. 

(a) Kick type A. It is a single-Maxwellian (for the sake of simplicity) distrib ution f(v) 



v^JjJp) exp(— (v /v-p) 2 ), as suggested by pulsar proper motion measurements (IHobbs et al. 



20051 ). The parameter v p is varied within the range 100-500 km s _1 . For NSs from single 
progenitors we always assume this type of the kick. 

(b) Kick type B. It assumes kick type A for NSs from single progenitors and from binary 
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Table 1. Model of the kick 



Model 


NS progenitors 


v p km s 1 


Initial mass 


SA 


single 


100 ... 500 


> 10M Q 


BA 


binary 


100 ... 500 


> 1OM 


BB 


binary 


30 
100.. .500 


8M -11M Q 
> 11M© 


BAS 

(BA + SA) 


50% binary 

+ 

50% single 


100. ..500 
100.. .500 


> 10M Q 

> 10Af Q 


BBS 


50% binary 


30 
100... 500 


8M© -11M© 
> 11M Q 



(BB + SA) + 

50% single 100... 500 > 10M Q 



progenitors with maximum masses during the lifetime above 11M . For NSs from binary pro- 
genitors with the maximum lifetime masses in the range 8-1 1M Q that e xperienced the Roche 



lobe overflow, in which the elect ron-capture core collapse can occur (IPodsiadlowski et al. 



2004 



van den Heuvel 



2004 



200^), v p is reduced and kept fixed at 30 km s 1 . 



According to the type of the kick assumed, we made calculations for different model 
populations of NSs, as summarized in Table [TJ Model SA means kick type A from single 
progenitors, BA means kick type A from binary progenitors, BB means kick type B from 
binary progenitors, BAS means that calculations are performed for equal number of single 
and binary progenitors with kick type A assumed, and BBS means that calculations are 
performed for equal number of single progenitors with kick type A and binary progenitors 
with kick type B. 

The kick type B appears to be more realistic and most of the results below are given for 
this model of the kick. It is necessary to note, that som e fraction of NSs origin ated from 



isolated progenitor can be formed after e -capture SN. 



Poelarends et al 



fbOQSh estimate 



this fraction as 4%, howeve r it can be sligh t 



is not well understood, yet. 



y hig her, because stellar physics prior to SN 



van den Heuvell (120071 ) notes that the fraction of normal radio 



pulsars originated from e -capture SNae cannot be high, because there is no room for large 
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fraction of low-velocity NSs among pulsars. On the other hand, one can speculaes that 
such NSs do not manifest themslves as normal pulsars. They can "hide" among radio silent 
compact objects. For example, close-by cooling NSs (aka "The Magnificent Seven") and 
central compact objects in supernova remnants. Pr oduct of e~-capture SNae are expected 



to be light. Light NSs cool slower (see a review in (jYakovlev &: Pethickl 12004)). And thos e 



Popov 



20081 ). 



types of radio silent NSs are know to be hot sources (see a brief review in 
Anyway, in this study we neglect the possibility that NS from initially isolated progenitors 
can be born after e~-capture SN. 

For both isolated and binary stars we also add some center-of-mass velocity, typical for 
massive stars in the galactic disc. For this velocity we use the Maxwellian distribution with 
v p = 10 km s _1 . The results of our calculations are found to be not very sensitive to the exact 



value of this parameter t o within a 
isolated massive stars see 



actor of ~ 2. 



Lequeuxi (1 1979( 1; IStond fjl99ll ). 



or a detailed study of space velocities of 



The rotational axes of both components are assumed to be aligned with the orbital 
angular momentum before the primary collapse forming the first NS in the system. The 
supernova explosion is treated in a standard way in terms of an instantaneous loss of mass 
of the exploding star. The effect of the kick on the post-explosion binary orbital parameters 



the description in e.g. 



Hills 


1983; 


Kalogera 


200o[ 



Grishchuk et al. 



20011 1 . The first supernova 



most likely occurs when the binary orbit is circular (unless the initial binary is very wide so 
that tidal circularization is ineffective), while the second explosion in the binary can happen 
before the orbit has been tidally circularized. Possible mass transfer phases before the second 
collapse (such as the common envelope stage and mass accretion onto the neutron star due 
to the Roche lobe overflow at the periastron passages) are assumed to effectively circularize 
the orbit. In the abs ence of mas s transfer the tidal evolution of the orbital eccentricity is 
treated according to IZahnl (119771 ) . When treating the SN explosion in an eccentric binary, 
the orbital position of the explosion is chosen randomly with a uniform distribution over the 
mean anomaly. 
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3 RESULTS: THE ROLE OF BINARIES 

In this section we discuss the effect of binary progenitors of NSs on the pulsar space velocities 
and spin-velocity misalignment. To this purpose, we calculate populations of isolated NSs 
arising from purely single and binary progenitors (models SA, BB), and populations resulted 
from equal number of single and binary progenitor stars (models BAS, BBS). 

3.1 Space velocity of young NSs 

First we consider the NS velocity distributions calculated from isolated and binary progen- 
itors. An isolated NS originating in a binary system acquires both the natal kick velocity 
and the orbital velocity of the exploding star immediately before the collapse. The result 
is sensitive to the mutual orientation of the kick and orbital velocities. In Fig. 1 we plot 
space velocities of NSs calculated for single and binary progenitors in comparison with a 
purely Maxwellian velocity distribution. As expected, the influence of the orbital velocity in 
binary progenitors is insignificant for large kick velocities (v p ~ 300 km s _1 , as was adopted 
in this plot). However, for kick type B, where a noticeable fraction of NSs in binaries re- 
cieved a reduced kick, the emerging NS velocity distribution slightly deviates from the pure 
Maxwellian law at low velocities. It appears quite difficult to check observationally such a 
small effect of e~-capture SN on the velocity distribution. 

3.2 Spin- velocity alignment in young NSs 

The effect of binary progenitors is more interesting for the expected NS spin- velocity align- 
ment. In Fig. [2] we plot the calculated spin- velocity alignment angle 6 for different kick-spin 
angles 6. The results are shown for young NSs in the model BBS. The spin and velocity 
vectors appear more aligned for a strong initial kick-spin alignment (i.e. small angles 8). 

The effect of binaries can be also illustrated (see Fig. [3]) showing the relative contribution 
of single and binary progenitors to the pulsars with given v and 5. Calculations are made for 
the models SA and BB with v p = 300 km s" 1 , and a fiducial value of the kick-spin alignment 
of 9 = 8° (the justification is given in Section 4). It is seen that NSs from single progenitors 
(model SA) are more abundant (i.e. the ratio of single to binary progenitors is more than 
one) at low velocities v (except for very low values < 30 km s _1 ) and small angles 5. In 
contrast, pulsars with v > 100 km s _1 and 5 > 15° - 20° are more likely to be born in a 
binary than from a single star (the single/binary progenitor ratio is smaller than one). 



Binaries and NS velocities 9 



0.0015 



0.0010- 



0.0005 



0.0000 



«— Model SA 
a Model BB 

Model BBS 
• — Maxwell 




50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 



V (km/s) 



Figure 1. The space velocity distribution of isolated young NSs (pulsars). The results are shown for the model SA (black 
squares), model BB (black triangles) and model BBS (grey upside-down triangles), with v p = 300 km s _1 . A pure Maxwcllian 
distribution (filled circles) is plotted for comparison. Note a small deviation from the pure Maxwellian distribution at low 
velocities for models BB and BBS. 




Figure 2. The spin-velocity alignment effect (the distribution of angles 8 between the spin axis and space velocity of a pulsar) 
for different kick-spin angles 8. Calculations are made for the model BBS with v p = 300 km s — , dN/dq ~ q 2 . The isolines on 
the plot stand for fO equal intervals of the distribution dN/dO/d8 from to the maximum value 2.4 ■ f0~ 3 . 
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Figure 3. The ratio of single and binary progenitors in the population of isolated NSs with different <5 and v. The calculations 
are shown in terms of N(model SA)/N(model BB) with v p = 300 km s — 1 and a narrow kick-spin angle 9 = 8°, dN/dq ~ q 2 . 
The values > 1 correspond to domination of single progenitors. It is seen that NSs from isolated progenitors are more abundant 
at small velocities v (except very low values < 30 km s — 1 ) and angles S. In contrast, pulsars with v > 100 km s _1 and 
5 > 15° — 20° are more likely to be born in a binary than from a single progenitor (the single/binary progenitor ratio is smaller 
than one). 



Let us consider more closely the spin-velocity alignment in pulsars originated from dif- 
ferent progenitors (Fig. Hj). A pulsar can result either from a single massive star or from the 
primary or secondary component of a massive binary system. In pulsars from single stars 
the spin-velocity alignment should reflect the kick-spin correlation during the NS formation. 
In pulsars from binary progenitors the initial correlation can be different depending on the 
conditions and history of the NS formation. The NS formed after the collapse of the pri- 
mary component Mi (NS1) can become an isolated radio pulsar (Mi — > NS1 = PSR1 I) 
immediately after the first supernova explosion if the binary becomes unbound, or after 
the second supernova explosion (Mi — > NS1 = PSR1 II) (filled squares and crosses in 
Fig. HI respectively). Pulsars can also result from the collapse of the secondary component 
M 2 when the binary system has already been disrupted by the first supernova explosion 
(M 2 — > NS2 = PSR2 I) or when the system became unbound after the second supernova 
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(M 2 — > A^S"2 = PSR2 II) . These cases are marked in Fig. [4] with filled circles and triangles, 
respectively. If the system remained bound after bo th explosions, a binary NS is formed; we 
discussed their coalescences in our previous paper (jPostnov fc Kuranovil2008l ). For compar- 
ison, open circles in Fig. H] show the spin-velocity angle distribution for pulsars originated 
from single progenitors. 

Fig. H] demonstrates that the population of single pulsars originated from binaries is 
largely dominated by those born after the disruption of the parent binary system by the 
first supernova (PSR1 I and PSR2 I), which constitute roughly ~ 50% and ~ 40% of the 
entire population, respectively. 

The effect of the e~-capture SNs, which is taken into acount in the models BB and BBS, 
should be mentioned explicitly. PSR1 I in these models are NSs originated from stars with 
Mi > 11 M , i.e. those which did not experienced an e~-capture SN. Conditions for an 
e~-capture SN are such that a system always survives after the first explosion. PSR2 II 
are mostly NSs from systems where the first explosion was an e~-capture event presumably 
producing small kicks so that the binary system survives, while the second explosion was due 
to the usual core collapse event. This happens when the secondary companion significantly 
increased its mass above 11 M Q due to mass transfer from the primary companion. 

Pulsars originated from binaries disrupted after the second supernova (PSR1 II and 
PSR2 II) provide a minor contribution but have, on average, larger spin-velocity alignment 
angles 5. Note that recycled pulsars, which can be spun-up by accretion in the binary after 
the first SN explosion, are among our population PSR1 II and constitutes less than few per 
cent of all pulsars. 
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Figure 4. The cumulative distribution of the spin- velocity alignment angles <5 for different subtypes of pulsar progenitors. From 
bottom to top: pulsars from primary components in binaries disrupted after the second explosion (PSRl II, crosses); pulsars 
from secondary components in binaries disrupted after the second explosion (PSR2 II, triangles); pulsars from secondary 
components in binaries disrupted after the first explosion (PSR2 I, filled circles); pulsars from primary components in binaries 
disrupted after the first SN explosion (PSRl I, squares). Stars show the total population of pulsars from all binary progenitors. 
Pulsars from single progenitors are shown by open circles. The calculations are shown for model SA and model BB with 
u p = 300 km s- 1 and 6 = 8°, dN/dq ~ q 2 . 



4 DISCUSSION 

It is interesting to compare the calculated spin-velocity angles 6 in single pulsars with ob- 



servations assumi ng different model pa ramet ers. Un 



measurements. In 



Ng & Romanil (120071 ) and 



Rankin 



ortunately, there are only a few precise 
(120071 ) the authors present lists of pul- 



sars for which the angle 5 can be inferred from observations. The cumulative (normalized) 
distributions of spin-velocity alignment angles (folded around 45° due to the 90° ambiguity 
in alignment angles measured by radio polariz ation) in pulsars fro m these lists are shown 
in Fig. [51 For several well-studied objects (see iNg &: Romanil 120071 and references therein) 
this angle can be determined from a detailed 3D-modeling with account of the X-ray data 
about the pulsar wind nebula shape, etc. For most pulsars the existing estimates of the spin- 
velocity angle 5 are based only on radio observations. In the last case the angle between the 
pulsar spin axis and its space velocity is derived from measurements of the proper motions 
of the pulsar on the sky and the polarization angle of radio emission. The angle determined 
as a difference between the proper motion direction and the pulse polarization angle spans 
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Figure 5. The cumulative distribution of the observed spin-velocity angles in pulsars folded around 45° due to the 90° 
ambiguity in alignment angles measured by radio polarization (data from Rankin 2007, filled squares, and Ng & Romani 2007, 
filled circles). 



the range from ~ —60° to ~ 100° ( IJohnston et al. 



2007 



Rankin 



20071 ). However, as we are 



interested just in the angle between the velocity and rotation axis, irrespective of the veloc- 
ity direction, this angle should lie bet ween and 90 degree s. In addition, projection effects 



complicate the interpretation of data (IJohnston et al 



well-studied objects, for example Vela (jNg fc Romani 



2005 



. On the other hand, for several 



20071 ). it is known that although the 



difference between the proper motion direction and polarization angle points to the spin 
axis perpendicular to the velocity, in reality 5 is close to zero. This is interpreted as be- 
ing due to the pulsar predominantly emitting ra diation which is polarized perpendicular to 



the magnetic field lines. So it is quite common (jNg &: Romani 



2007 



Johnston et al. 



2003) 



to fold data in such a way that the value of the spin-velocity angle lies within the range 
0° < 5 < 45°. 

The spin- velocity angle distributions are calculated for different kick-spin angles 9, differ- 
ent kick velocities and assumptions abo ut the initial binary mass ratio dist ribution. The com- 



parison is made with both data sets by (jNg fc Romani 



2007 



Rankin 



20071 ) . The Kolmogorov- 



Smirnov (K-S) test was applied to test the nul hypothesis that the observed (Fig. [5]) and 
calculated alignment angles are realizations of one distribution. The results are shown in 
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Figs. [6] [HI The calculated and observed spin- velocity angles S were binned in 23 two-degree 
intervals to increase the number of bins above 20, as recommended for the K-S test. Since the 
observational data is folded back to 5 < 45° the same folding was applied to the calculated 
alignment angles. 

Fig. [6] shows a 3D graph for the K-S test of the model distribution as a function of 
the fraction of binary progenitors (we mix models SA and BB) and the kick-spin angle 9. 
Both panels of Fig. [6] are shown for v p = 300 km s -1 and dNjdq ~ q 2 . Th e left and righ t 
panels present comparison with data from Ng fe Romani ( 2007 ) and from Rankin (2003), 
respectively. It is visible that better fits are obtained for a high fr a ction of binary systems. 
Especially, in the case of comparison with the data from iRankinl (120071 ). small fraction of 
binary progenitors is disfavoured. 

We studied the influence of the fraction of binary progenitors for different kick velocities 
on the quality of the K-S test. For both samples of observed pulsars (from Ng & Romani 
and from Rankin) for reasonably high velocities (v p > 100 km s _1 ) the fraction of binary 
progenitors should be higher than ~50% to produce high quality of the test. For more 
realistic velocities (v p ~ 200 — 300 km s _1 ) comparison of our calculations with observational 
data favours a high fraction of binary proge nitors >70%, in correspondence with observations 
of binary frequencies among massive stars (IMason et al.lll998 ; 



Gies 



2008|). 



Fig. [7| shows the K-S test calculated for model SA (upp e r pan els) and BB (lower 



as a 



unction of angle 9, for pulsars from 



panels) , 



Ng fc Romani! (120071 ) (to the left) and 



Rankin 



(120071 ) (to the right). Fig. [8] shows the K-S test for models BAS and BBS with different 
assumptions about the initial mass ratio distribution f(q). 

It is seen from Fig. [7] -[8] that the K-S test in all cases does not show good correspondence 
between the observed and calculated pulsar spin- velocity distributions, which may be related 
to a poor statistics of pulsars with known 5. Nevertheless, the K-S maxima are clearly 
distinguished in all plots, so the KS-test can be used as a guide to assess the relative goodness 
of the fits. 



For all kick models and initial mass ratio distributions, the data taken from 



Ng fc Romani 



( 120071 ) are best-fitted by kick velocities co nfined within the angle 9 ~ 20°. In contrast, the 
comparison with a more numerous data by IRankinl (120071 ) favours a narrower natal kick-spin 
alignment with 9 ~ 5 — 10°. This ex plains our cho ice of the fiducial value 9 = 8° we used 
in Figs. [31 HI Note also that data by IRankinl (120071 ) can not be fitted by single progenitors 
only (model SA, right upper panel in Fig. [7]). 
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Figure 6. The Kolmogorov-Smirnov test for the observed and calculated spin-velocity alignment angles 8 in pulsars as a 
function of the assumed kick-spin alignment angle 9 and the ratio of binary and single progenitors of the population of isolated 
NSs. The latter is given as N(model BB)/[N(model BB)+N(model SA)]. Calculations are shown for the model BBS. The initial 
binary mass ratio distribution is f(q) ~ q 2 , and the kick velocity amplitude is equal to 300 km s — 1 . Left panel: comparison 
with the data from Ng & Romani (2007). Right panel: comparison with the data from Rankin (2007). 

Fig. [8] clearly shows that the results only slightly depend on the chosen kick model and 
the kick velocity amplitude, and are virtually independent of the assumed initial binary mass 
ratio distribution f(q). 

It may seem that the actual shape of the kick velocity distribution can be an important 
parameter. In our study we restrict ourselves to considering only the simplest case of a 
single-peak Maxwellian distribution for most of NSs is in agreement with observations. As 
our kick model B assumes the formation of a certain fraction of low-kick NSs and, hence, 
low-velocity pulsars, the obtained insignificant dependence of the calculated spin-velocity 
alignment on kick models (at least with the current pulsar data at hands) suggests the true 
shape of the kick being not so important. However, it does matter for the estimation of 
the relative fraction of single and binary pulsar progenitors at given 5 and v (Fig. [3]). For 
example, if the true kicks can take values close to zero even for single progenitors, it would 
be hard to distinguish binary and single progenitors from such an analysis. 
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Figure 7. The Kolmogorov-Smirnov test for the observed and calculated spin-velocity alignment angles 8 in pulsars as a 
function of the assumed kick-spin alignment angle 9. Calculations are shown for the initial binary mass ratio distribution 
f(q) ~ q 2 and models SA and BB. Squares, circles and triangles mark the kick velocity amplitudes 100, 300 and 500 km s — 1 , 
respectively. Left panels: comparison with the data from Ng & Romani (2007). Right panels: comparison with the data from 
Rankin (2007). 



5 CONCLUSIONS 

In the present paper we considered the role of binary progenitors of neutron stars in the 
apparent distribution of space velocities and spin-velocity alignment of young pulsars. We 
used a Monte-Carlo synthesis of pulsar population from single and binary stars under dif- 
ferent assumptions about the NS natal kick model (form, amplitude and possible reduction 
of the kick in binary porgenitors with masses from the range 8-11 M ) and initial binary 
mass ratio distributions. We compared the calculated spin-velocity alignment distributions 
with observational data obtained from radio polarization measurements. We arrived at the 
following conclusions: 

1) . The observed space velocity of pulsars is mostly shaped by the assumed natal kick 
velocity form and amplitude (Fig. [T]). The effect of binary progenitors is negligible. 

2) . The possible kick-spin alignment during the formation of a NS strongly affects the 
spin-velocity correlation in pulsars. For the kick model considered, single progenitors are 
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Figure 8. The Kolmogorov-Smirnov test for the observed and calculated spin-velocity alignment angles S in pulsars as a 
function of the assumed kick-spin alignment angle 9. Calculations are shown for the initial binary mass ratio distribution 
/(?) ~ Q° an d /(<?) ~ Q 2 (the upper and lower panels, respectively) and different models (BAS, BBS). Squares, circles and 
triangles mark the kick velocity amplitudes 100, 300 and 500 km respectively. Left panels: comparison with pulsar data 

from Ng & Romani (2007). Right panels: comparison with pulsar data from Rankin (2007). It is seen that the form of the 
initial binary mass ratio distribution and type of kick velocity models (BAS, BBS) have no influence on the results. 
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more probable for pulsars with space velocities > 50 km s _1 and tight spin-velocity align- 
ment (< 10°) (Fig. [3]). Binary progenitors are favoured for pulsars with 5 > 15° and 
v > 100 km s , and for low velocity pulsars (< 30 km s _1 ) for all values of 5. 



3) . The compari son with current measurements of pulsar spin- velocity angles (INg fc Romani 



2007 



Rankin 



20071 ) does not allow us to distinguish between different kick models - the 



agreement between the calculated spin-velocity alignment and the observed distributions is 
not very good (Fig. [71 [H])- Nevertheless, the assumption of natal kick-spin alignment during 
NS formation appears to be more important than the kick veloc ity distr i butio n model. 
Single progenitor s only (model SA, Fig . [7J do not fit data by iRankinl (120071 ) and provide 



worse fit for data by 



Ng fe Romanil (120071 ) than binary p rogenitors. Bin a ry pro genit ors only 



Rankin 



(mode l BB, Fig. [7J can equally well fit both data by INg fc Romanil (120071 ) and 
( 120071 ) for the kick-spin velocity angle 9 ~ 5 — 20°. Models with equal numbers of single and 
binary progenitors of isolated young NSs (models BAS, BBS) equally well fit observations 
irrespective of the assumed form of the initial binary mass ratio distribution f(q) (Fig. [8]). 
Models without binary progenitors fit observations significantly worse for standard mean 
kick velocity amplitude v v ~ 300 km s . 
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